Transport in molecular states language: 
Generalized quantum master equation approach 
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A simple scheme capable of treating transport in molecular junctions in the language of many-body 
states is presented. An ansatz in Liouville space similar to generalized Kadanoff-Baym approxima- 
tion is introduced in order to reduce exact equation-of-motion for Hubbard operator to quantum 
master equation (QME)-like expression. A dressing with effective Liouville space propagation similar 
to standard diagrammatic dressing approach is proposed. The scheme is compared to standard QME 
approach, and its applicability to transport calculations is discussed within numerical examples. 

PACS numbers: 85.65.-|-h 85.35.Be 73.63.Kv 73.23.Hk 



I. INTRODUCTION 

Quantum transport in nanoscale systems is on the fore- 
front of research in many fields of study. In particu- 
lar, progress of experimental capabilities in the field of 
molecular electronics brings new theoretical challengesii 
Resonant transport with strong on-the-bridge interac- 
tions is one of them. It is probably the most impor- 
tant regime for possible future applications (e.g. for logic 
and memory molecular devices). Unlike usual mesoscopic 
systems, molecular electronic (and vibrational) structure 
may be very sensitive to reduction/oxidation. Thus res- 
onant transport in molecular junctions has to be de- 
scribed in the language of molecular states (states of the 
isolated molecule as a basis for description of reduced 
non-equilibrium dynamics), rather than in the language 
of (effective) single-particle orbitals, which is generally 
accepted in the molecular electronics community. Re- 
cent experiments on simultaneous measurements of cur- 
rent and optical response of molecular junctions^ make 
need for such formulation even more pronounced, since 
molecular states is a natural language of all the (equilib- 
rium) molecular spectroscopy. Besides, molecular states 
based formulation of transport makes it potentially pos- 
sible to incorporate standard quantum chemistry molec- 
ular structure simulations as an input to transport cal- 
culations. 

Necessity of many-body states type description of 
transport in molecular systems has been realized. 
Among the approaches one can mention scattering 



theory^ (generalized) master equation 



4.5.6.7.8.9 



density 



matrix j-'^^'-'^-'^i-'^^i-'^'^i-'^^'-'^^i-'^^ and non-equilibrium Green func- 
tion (NEGF) based schemes i^^i^°'^^ Each of these 
schemes has its own limitations. Scattering theory when 
applied to transport problems disregards junction char- 
acter of molecular system which may lead to erroneous 
predictions f^ii^ Standard formulations also miss crucial 
physics, e.g. effective attractive electron-electron inter- 
action via phonons (bipolaron formation), energy ex- 
change (heating and cooling effects) between successive 
tunneling events, target distortion due to quasi-bound 
states etc. Master equation (or generalized master equa- 



tion) approaches are used quite often to describe hop- 
ping transport, i.e. situation when correlations in the 
system (both in space and time) die much quicker than 
electron transfer time (fixed by contact/molecule cou- 
pling). Besides, they become inadequate in the off- 
resonant tunneling (super-exchange) situation. Density 
matrix schemes usually formulated within quantum mas- 
ter equation form, often miss broadening of molecular 
states due to coupling to the contacts and coherences re- 
sponsible, e.g., for elastic channel renormalization at the 
inelastic threshold (see also discussion below). NEGF 
based approaches in the language of molecular states 
are among the most advanced methods for treating non- 
equilibrium molecular systemsi^ Two important draw- 
backs of the approach are its complicated character and 
absence of proper commutation relations for the Hub- 
bard operators. The first means that applicability of the 
method is limited to simple cases only. The second may 
lead to unphysical consequences (e.g. non-Hermiticity 
of the reduced density matrix) at approximate level of 
treatment i^>^ 

The goal of the present paper is to formulate an ap- 
proximate scheme for treatment of transport in molec- 
ular junctions in the language of molecular states, ex- 
ploring connection between Green function and den- 
sity matrix based approaches to transport. We start 
from the NEGF-like consideration and derive QME-like 
equation, pointing out approximations involved in the 
derivation. Note that similar approach within a single- 
particle orbitals language was used in Ref. [H. Note also 
that in our consideration we go beyond strictly Markov 
limit of Ref. [l^. We work with Hubbard operator as 
a natural object capable of describing excitation in the 
molecule as transitions (we restrict our consideration to 
single-electron transitions) between many-body molecu- 
lar states. First we show how exact equation can be re- 
duced to QME by introducing ansatz in Liouville space 
similar to the generalized Kadanoff-Baym approximation 
(GKBA)^ in Hilbert space. Second, we identify dia- 
grams on the Keldysh contour corresponding to the pro- 
cesses described by QME, and in the spirit of Green 
function diagrammatic techniques dress this diagrams 



2 



by effective Liouville space dynamics. Tfie last is ob- 
tained from the exact equation within Markov approxi- 
mation. Section introduces molecular junction model 
and presents equation-of-motion for the Hubbard oper- 
ator. In section IIIII we introduce an ansatz in the Li- 
ouville space, which being similar to the GKBA, allows 
reduction of the exact EOM for the Hubbard operator to 
QME. Here we also discuss a dressing procedure. Sec- 
tion IIVI presents analytical consideration for the simple 
resonant level model, and discuss general numerical pro- 
cedure. Several numerical examples are presented as well. 
Section W\ concludes. 



II. MODEL 

We consider molecular junction which consists of two 
contacts (L and R) coupled through the molecule (M). 
Contacts are assumed to be reservoirs of free electrons 
each at its own equilibrium. All the non-equilibrium 
physics takes place at the molecule. Hamiltonian of the 
system is 



H = Ht 



H 



M 



Hr 



V = Hq + V 



(1) 



where Hm is a full molecular Hamiltonian, i.e. Hamilto- 
nian of isolated molecule with all on-thc-molecule inter- 
actions included, V is molecule-contacts coupling 



V = 



E 

m£ALke{L,R:} 



(2) 



and Hk [K = L, R) represents contacts 



(3) 



keK 



Here (cfc) and dl^ (dm) are creation (annihilation) op- 
erators for electron in a state k in contact K and state 
m of molecular Hamiltonian, respectively. 

We introduce many-body states of isolated molecule 
\N,i > with N being number of electrons on the molecule 
and i standing for a set of all other quantum numbers 
characterizing particular state of the molecule in the 
charging block N. These states are assumed to be or- 
thonormal 



< N,j\N',i' >= S_ 



N,N' Oi.i' 



(4) 



Note, that generalization to nonorthogonal basis is avail- 
able in the literatures^ but we will stick with orthonor- 
mal basis in order to keep notation as simple as possible. 
Molecular transitions (in our case due to coupling to the 
contacts) are naturally described in the language of Hub- 
bard operators 



X(N,^,N',^') = \N,i><N',l'\ 



(5) 



One of important (in our case) transitions is oxida- 
tion/reduction of the molecule by one electron, i.e. tran- 
sition between neighboring charge blocks 

M = {N,i]N + l,j) M = {N + l,j]N,i) (6) 

In terms of these states molecular Hamiltonian is 



N,i,3 



N,i,3 



(7) 

If the many-body states are chosen to be eigenstates 
of the molecular Hamiltonian, H^^^ = E^'^^Sij with 

) being energy of the molecular eigenstate \N,i >. 
Molecule-contacts couping, Eq.®, becomes 

V=Y1 {VkMclXM + VMkXM^k) (8) 



k,M 



where 



VkM = J2 ^/^^ < ^' + 1, j > 



(9) 



t 

Al- 



and V^^, = V^j^. Note, that = X 

In our previous publicationj^i we considered applica- 
tion of a method originally formulated in Ref. to in- 
elastic transport in molecular junctions. The main object 
of interest in this consideration was many-body (Hub- 
bard) Green function on the Keldysh contour 

G(a;6),(c;d)(r,r') = -i< T,Xab[T)Xl^{T') > (10) 

where a, b, c, d are many-body states of an isolated 
molecule, Tc is the contour ordering operator, and r, 
t' are contour variables. The consideration leads to a 
formulation similar to standard diagrammatic technique, 
with series of functional derivatives in place of expansion 
in small parameter for conventional diagrammatic con- 
sideration (for a detailed discussion see Ref. [2l|). Thus 
obtained machinery is quite general, but may be too 
heavy for realistic applications. Also at approximate level 
of treatment it may lead to non-physical results.— 

Here the main object of interest will be the opera- 
tor Xab{t) = e^^^Xabe'""^*^, where \a >= \Na,Sa > and 
|6 >= \Nb,Sb > are many-body states defined in dH), and 
t is time. Our goal is utilizing Green function techniques 
find a (approximate) connection to density matrix based 
considerations (in a manner similar to that of Ref. [27l ). 
and use resulting scheme as a simplified version of a pro- 
cedure considered in e.g. Ref.[2ll. Note, that we are going 
to go beyond standard QME considerations of transport 
(see discussion below). Our starting point is EOM 



< Xabjt) > 

dt 



H; Xab{t) 



(11) 



Taking commutator in the right side of pT|) yields cor- 
relation functions of the form (for detailed derivation see 
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Appendix El < )(t)cfc(i) > and < 4(i) -^(...)(*) >• 
As usua l^^'^^ these correlation functions can be treated 
as lesser projections of Green functions 

GcxiT,T') = -^ < r,Cfc(T)xt )(t') > (12) 

GxciT,T') = -^ < r,l(...)(T)4(T') > (13) 

respectively, taken at equal time. The last can be ob- 
tained by applying Langrcth projection rules^ to on-the- 
contour EOMs for ^ and ^ 

GexiT,T')^J2 I d^i9kiT,n)VkMGM,..Xri,T') (14) 
Gxc{t,t') = Y, I dTiG...M{T.Ti)VMk9k{TuT') (15) 



M 



where Gm,... {G...,m) is defined in ([rU]) . VkM {VMk) is 
introduced in ([9]), and 



gk{T,T') = -t <T,Ck{T) 4(t')> 



(16) 



is Green function of free electrons in the contacts. 

Using lesser projections taken at equal times of p4] 
and (|15|) in (fTTj) leads to (see Appendix [Al for details) 



d < Xgbjt) > 
dt 



< X 



{Na,s-Nb,si, 

t 



{t)> 



{ ^(Ar„,s„;Afb + l,s),Al(*'*l) ^A1,(Afb,Sb;Ar6 + l,s)(*l ~ 



(-1) 



(17) 



~ ^riVa-l,s;iV6,Si,),A4(^'*l) ^A1,(Wa-l,s;7V„,s„)(*l ~ 



Here I]^^^^(t) are greater and lesser molecular self- 
energies due to coupling to the contacts 

^MlM.^t)- E ^MlXit) (18) 
K=L.R 



(19) 



with 5^ ' (i) being greater and lesser projections of ([16] 



g>(i)=-z[l-n,.]e— 



(20) 
(21) 



and G^^^, (ti,t2) are greater and lesser projections of 

m 

Gr.;fc)^(.;,)(ii,i2) < XMXl,{h) > (22) 

Gf„;fc),(c;d)(*i-*2) =±^ < ^t,(t2)X,fc(ti) > (23) 



Note, in (|23p ' + ' stands when both A^i and A^2 are 
transitions of Fermi type, and ' — ' otherwise. For future 
reference we also define a damping matrix in Liouville 
space 



^ Mi,M2 — 



(24) 



Expression for the current can be derived in a similar 
way (see Eq.(lO) of Ref. [2l|) 

= l E r dt' 

+ G>^^^,{t,t')^<^,^^{t' -t) (25) 

^ ^A1,A1'(^ ^ ^ ) G%^i _x^{t , t) 

Equations ([T7]) and are exact, however their right 
sides are expressed in terms of Green functions. Our goal 
now is to introduce approximate scheme in order to close 
(fT7)) in terms of < Xab{t) >, thus finding connection to 
QME. This approximation is introduced and discussed in 
the next Section. 



III. GENERALIZED QME 

Before introduction of an ansatz to close Eq. p7|) we 
note close connection between EOM for Xab{t) and den- 
sity matrix element Pha{t)- Indeed, 



Pba{t) =< X6a|e~*"|po X6a|e-'"/5o > (26) 

= < e'^''^Xba\P() » = < Xab{t) > 

where C is the total LiouviUian and < A\B >= Tr[At B] 
is scalar product in Liouville space. Hence we expect that 
Eq. (|17p after introducing approximation expressing it in 
terms of Xab only should result in QME. 

Correlation function of the type (|22p can be exactly 
written in Liouville space as 

< Xabitl) Xl^{t2) >= (27) 

e{ti - t2) « XbjK\c-"^^''-'''^\XdcP{t2) > 
+ e{t2~ti) « X,jK\e-'^^*'-*'\p{ti)Xab » . 



We introduce the projector superoperator 



(28) 
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which disregards nonequiHbrium features in the leads and 
decouples system and bath dynamics. The ansatz, we 
propose, replaces (pTj) by 



< Xabih) xUh) >^ (29) 

0{h - t2) « XbjK\c-'^'''''''MXdcP{t2) > 

+ e{h - h) < x,jK\c-'^'-*'-''^v\p{ti)Xab » . 

Next we introduce retarded and advanced Green func- 
tions in the Liouville space (see Appendix IB|) 

G:j,ranit) = -*^(i) « X,./,^ je"*^* |X„„,p^' > (30) 

= -ie{t) < Xji\Ueff{t)\X„„. > 

G:j,ran{t) ^ ^S{~t) « X^J^ » (31) 

= t9i-t) < X,,\Ulfji-t)\Xn„^ > 

where the effective propagator in the molecule space 
reads 



Using (|30p and (PT|) . we can rewrite as 

<l,fc(il)lt^(t2) > = 

[^a6./e(il - ^2) < X/,(t2)l,yt2) > 

~<X,fc(ti)l]^(tl)>ae%d(tl-i2) 

= «E [^a6,rnd(il - ^2) < X^c{t2) > 
m 

- < Xamih) > Qmhcdi^l ^ ^2) 



(32) 



(33) 



where second equality follows from orthonormality con- 
dition (g]). 



Similar consideration for correlation function (j23p 
leads to 



<Xl{t2)Xab{h)> = 
*E [^abjeitl - h) < XUh)Xfe{t2) > 

- < Xlj.{h)Xabitl) > g^f,cditl~t2) 

= *E [S:b,cmih-t2)<Xd,nit2) > 
rn 

- < X„ib{tl) > G^m.cdi^-l - *2) 



(34) 



It is interesting to note that (|33p and ([34)) can be con- 
sidered as the Liouville space analog of the generalized 
Kadanoff-Baym ansatz.— 



Using (|33p and ([34]) in (fTTl) closes the latter equation 
in terms of DM pba{t) =< Xab{t) > only 



dpujt) 
dt 



^E 1 ^Ni,N35N2,Ni E ('^'2: 
3,4 I s 



E 



4-oc 



^(2;JVi + l,si)(4;Ar3 + l,S2)(^ ^1 )^(3;Ar3 + l:S2)(l;Ari + l,si) (^1 



^(2;Ar2 + l,si)(4;Ar4 



+ 1,S2)(^ ~ *l)^r3;Af4 + l:S2)(l;JV2 + l,si)(*l - 1) 
~ ^(JV2-l,si;l)(Af4-l,S2;3)(* ~ *l)^^4-l,S2;4)(7V2-l,si;2)(^l ~ 
+ ^(Ari-l,si;l)(7V3-l,s2;3)(* ^ *l)^{V3-l,s2;4)(Ni-l,si,2) 

(tl - i) 

(^1)^' ^ (^('w2-l,si;l)(4;Af3 + l,S2)(^ ~ *l)^(3;W3+l,S2)(Af2-l,si;2)(^l ~ 

^ ^(}Vi-l,si;l)(4;JV4 + l,S2)(^ ^ ^l)^(3;JV4 + l,S2)(Afi-l,si;2)(*l ^ 1) 
~ !='(2;7Vi + l,si)(JV4-l,S2;3) l,S2;4)(l;Wi + l,si) 

(il - t) 



^(2;A'2 + l,Si)(iV3- 



-1,S2;3)(* *l)^('V3-l,S2;4)(l:Af3 + l,si) (ti-i))] Xp34(il)} 



(35) 



This is a generalized non-Markovian QME. Note, that prcfactor (— 1) 



Ni-N-2 



coming from coherences between 
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different charge blocks is usually lost in the standard 
QME derivations. 



To make ((35)) more tractable below we assume 
Markovian generator £e// (e.g. Markovian Redfield 
generatori2ii^i2i) for retarded and advanced Green func- 
tions §UI and §11 



Ueffit) « e 



(36) 



So that 



-iCe f ft 



(37) 



Gr^.mnit) ^ ^^{^t) « X,. | e"'''^^'^* |^„™ » ■ (38) 



The ansatz (PH]) together with ([571) ^i^^ (EH) is equivalent 
to use of the regression formula on the Hubbard Green 
function. This procedure is commonly used to calculate 
multipoint correlation functions using effective Marko- 
vian propagators^ii^i^^. 

The standard non-Markovian QMEji^ is obtained from 
([55)) by using in ([57)) and ([55)) the free molecular evolution 
£m = [Hm,-] instead ofthe effective one £e//. Note that 
difference between standard and generalized versions of 
QME is similar to dressing of diagrams in GF diagram- 
matic technique. Note also that the standard QME by 
itself can not reproduce, e.g., broadening of molecular 
levels due to coupling to the contacts as noted in Ref . [l^. 

Below we use the Markovian Redfield equation to get 
jCe// (see Appendix [C)) . Its spectral decomposition 



^eff - ^ |i?7 > < 



(39) 



with eigenvalues and left \L^ 3> and right \R^ 3> 
eigenvectors, provides a numerically tractable scheme to 
deal with generalized QME ([55)) by utilizing 
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(40) 

7 

(41) 



Steady-state for ([55)) is given by the right eigenvector 
with zero eigenvalue of the Liouvillian corresponding to 
the Markov hmit of ([35l). 



of < X( ) > can be obtained using (|33p and ([34)) in ([25)) 



E 



2Re S('Ari,ii;Afi + lji),(A'2,i2:e)(^ ~ ^l) 

^ ^^2,i2;W2 + lj2),(^'i:»i;^i + lji)'-^l ~ 
X < l(jV2 + l,i2;e)(il) > (42) 
+ ^('7Vi,Ji;A'i + lji),(e;JV2-l-lj2)(* ~ ^l) 

X < -'^(e;Ar2,i2)(*l) > 



IV. RESULTS AND DISCUSSION 

As a first example we consider a simple resonant level 
model. One has two charge blocks (occupied and unoc- 
cupied level) with one state in each of them: |0 > and 
|1 >. The molecular Hamiltonian is Hm = |1 > eo < 1|- 
Current ([1^ in this case becomes 

lK{t) = jI' dh { (43) 

Ki,oi(^ - ti)Y.>{ti - t) - Y>{t - ii)go"i,oi(ii - t)] 
X Piiih) 

xpoo(ii)} 
where 

Goimit) = -^Oit)e-^^^"-^^/'^' ^ G^{t) (44) 
Qoimit) = ie{~t)e-'^'°+'^'^^' = C'it) (45) 

and r = Y.K=L.R^mm with defined in ([24]). 
Generalized QME ^} yields 



d/9ii(t) (i/9oo(i) 



/ dti { 



(46) 



dt dt 
[aSi.oi(i - ii)S>(ti - t) - S>(t - ti)eo"i.oi(ii - t)] 
X Pii(ii) 

+ Ki,oi(i - ii)S<(ti - i) - S<(t - ti)go"i,oi(ii - i)] 
xpoo(ii)} 
At steady-state ([46)) yields 

Pii = 1 - Poo = »^o (47) 
with no average occupation of the level 



^fL{E) + ^fR{E) 



Similarly, approximate expression for current in terms 



A{E) = 



(i?-eo)2 + (r/2) 



(48) 
(49) 
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where A{E) is spectral function and fniE) is Fermi dis- 
tribution in contact K = L,R. Using (gi]), (gSl), and (gT]) 
in ((43|) leads to the Landauer expression 



Ik = 



dE TLTf 



27r 



(50) 



Note, that generalized QME approach takes level broad- 
ening into account in a natural way contrary to the stan- 
dard QME considerations. 



U 



-a 

> 




0.0 0.5 1.0 1.5 
V (arb.u.) 



2.0 



FIG. 1: (Color online) Comparison of generalized (solid line, 
red) to standard (dashed line, blue) QME results for resonant 
level model, (a) Current vs. bias for single level, (b) Con- 
ductance vs. bias for single level coupled to a vibration. See 
text for parameters. 

Now we present several numerical examples. Figure [T] 
compares results of calculation within our generalized 
QME (solid line) and standard QME (dashed line) ap- 
proaches. Fig. [T^ shows current-voltage characteristic of 
single resonant level Eq model. Generalized QME ac- 
counts for level broadening due to coupling to the con- 
tacts, while standard QME approach misses the broaden- 
ing altogether. Parameters of the calculation are £o = 1, 
= = 0.1, Ef = m = 0, A^l =Ef + \e\Vsd. Here 
and below we use arbitrary units. Fig. [T]d shows conduc- 
tance vs. bias for the model of single level eg coupled to 
a vibration ljq. Once more, while generalized QME pro- 
vides reasonable results (compare eg to Fig. 4 of Ref.llSha 



standard QME approach is capable of prediciting only 
positions of the peaks. Parameters of the calculation are 
Tl = Tr = 0.05, cjo = 0.2, and M = 0.2. The last 
is strenth of electron-vibration coupling on the bridge 
with corresponding Hamiltonian M{a -f a^)fiQ, where a) 
(a) are creation (annihilation) operators of vibrational 
quanta and no operator of the evel population. Other pa- 
rameters are as in Fig. [IJi. Note, that in simulations we 
used small but finite broadening for the standard QME 
approach in order to avoid delta-function divergencies in 
conductunce. We also scaled the standard QME result 
in Fig. [T]d for convenience. 

Figures [2] present conductance maps for a quantum dot 
obtained within the generalized QME approach and sim- 
ilar to those obtained within many-body Green function 
technique (see Ref. [2]] ). Parameters of the calculation 
are level positions = —0.5 [a = {t,i}), molecule- 
contacts coupling Tk,c7 = 0.01 [K ~ {L,i?}), on-site 
repulsion [/ = 1, Fermi level Ep ~ 0. Electrochemical 
potentials in the contacts are fJ^L = Ep + |e|T4d/2 and 
/iij = Ep — \e\Vsd/'2'- Deviations from this 'standard' set 
for each calculation are specified below. Fig. [2^ shows 
conductance map for QD with level degeneracy removed 
by e.g. external magnetic field = —0.6 and = —0.4. 
For discussion on origin and intensity of peaks see Ref. Issl. 
Fig. [2)3 shows conductance map for QD coupled to a vi- 
bration loq = 0.1 and M = 0.1. In addition to elastic 
peaks vibrational sidebands corresponding to resonant 
inelastic tunneling are reproduced as well. Fig. [2j: shows 
conductance map for QD with asymmetric coupling to 
the contacts Tl^ = 0.01 and Vn^a =0.1. This result is 
similar to the one presented in Fig. 4 of Ref. [l^ 

Note, that vibrations in both Fig. [1)3 and[2)3 were intro- 
duced, as is usually done in resonant inelastic transport 
coniderations, with the help of small polaron transforma- 
tion. So that vibrational features in electron transport 
stem from the Franck-Condon factors calculated under 
assumption of unperturbed thermal distribution of vi- 
brational population. Actual vibrational states are not 
included in the current consideration, and their incorpo- 
ration into many-body state description will be described 
elsewhere. 

Finally, we consider a model of two-level bridge with 
coherences in the eigenbasis of the bridge induced by cou- 
pling to the contacts. This models was previously consid- 
ered in Ref . within standard QME approach. Figure[3] 
presents comparison between standard and generalized 
QME approaches. Parameters of the calculation are sim- 
ilar to those in Ref.— - eigenenergies of the bridge are 
£1 = 5eV are £2cV, strength of their coupling to con- 
tacts is = 0.3eV, = 0.2eV, = 0.4eV. 
For temperature we take physically reasonable value of 
T = O.OSeV. Figures [SK andjSb show current and one of 
the probabilities (probability of the system to be unoc- 
cupied) vs. applied bias. One sees that broadening due 
to coupling to the contacts is preserved in our scheme. 
Note, that broadening presented in Ref. [3 was due to 
unphysically high value of temperature chosen. Fig. [5J; 
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FIG. 2: (Color online) Conductance vs. applied bias Vsd and 
gate voltage Vg for a quantum dot (QD) within generalized 
QME approach. Shown are results for models of (a) QD with 
level degenraucy removed (e.g. by applied magnetic field), 
(b) QD coupled to a vibration, and (c) QD with asymmetric 
coupling to the contacts. See text for parameters. 



demonstrate influence of broadening on coherences (in 
local basis). Here we bring the two eigenenergies closer 
to each other, ei = 3eV, in order to make coherences 
due to coupling to the contacts more pronounced. One 
sees that taking level broadening into account changes 
the coherences essentially. 




^ 0.0 



o 
o 



-0.1 



2 3 
Vsd (eV) 



FIG. 3: (Color online) Two-level bridge with coherences in the 
eigenbasis of the system induced by coupling to contacts.— 
(a) Current and (b) probability for the system to be unoccu- 
pied vs. aplied bias - generalized (solid line, blue) and stan- 
dard (dashed line, red) QME considerations, (c) Real (solid, 
blue and dashed, red) and imaginary (dotted, blue and dash- 
dotted, red) parts of coherences in the local basis vs. applied 
bias for generalized and standard QME treatment, respec- 
tively. See text for parameters. 
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V. CONCLUSION 

Necessity for description of molecular transport in the 
language of many-body (isolated molecule) states, essen- 
tial for description of resonant tunneling and for study 
of optoelectronic devices, has been realized and several 
approaches were proposedi^'^'^i^i^ '^^i^^'^^i^'-'i^^ Here we 
introduce a simplified version of the Hubbard operator 
Green function approach considered in application to in- 
elastic transport in our previous publicationi^ The sim- 
plified approach is formulated for density matrix instead 
of GF and provides more easy way for calculating both 
timc-dcpcndcnt and steady state transport in molecular 
junctions. Starting from GF-typc consideration we intro- 
duce Liouville space analog of the generalized Kadanoff- 
Baym ansatz, which allows us to derive generalized QME. 
The latter differs from the standard QME by incorpo- 
rating effective propagation in place of free evolution. 
The procedure is similar in spirit to diagrams dressing in 
GF diagrammatic techniques. Capabilities of the scheme 
are demonstrated within model calculations. Application 
of the approach to opto-clcctronic response of molecular 
junctions is a goal for future research. 
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APPENDIX A: DERIVATION OF EQ. ([T7)) 

We start from Eq. (fTTI) , which after evaluating the com- 
mutator becomes 



d < Xabit) > 



dt 



< X 



(JV„,s„;Afb,s) 



(t)> 



{Na.s;Nt,st) 



(Vk,{Nt.Si-Nb + l,s) < 4(0^(^V„,.s„;JVi.+l,s)(^) > 



-V( 



V, 



k,{N^-l,s;N^,s. 



< cli*) ^{Na-l,s;N,,s,)it) > I 



where ... is sum over molecular states within charge 
block, ... is sum over states in the contacts, and fac- 



tor (—1)^" results from commuting Xab with Ck (c|). 

Correlation functions in the right of Ea. (|Al[) can be 
identified as lesser projections of the GFs (fT2|) and ((T3|) 
defined on the Keldysh contour. EOMs for this GFs are 
presented in and Taking lesser projection of 

the EOMs and applying the Langreth rules^ yields, e.g. 
for the first correlation function in (|Aip 

/ + 00 
dt, [(-1) 

{Na,Sa;Nb + l,s){t) > 5fc(^l - t) (A-2) 
+ e{t-h) (<X(jv,„,„;Ar, + i,,)(t) X|^(tl) > 

-(-1)^"'^^-' < ^Xl(^l)^(iV„...;^. + l..)(0 >) 

where g'^'^{t) are advanced and lesser projections of the 
GF (|16p . < ... >= Tr[. . . po] with initial density matrix 
taken as usual at infinite past, and where general prop- 
erty of GFs ^{1) = e{t)[G> (t) - G< {t)] was used for the 
Gxx GF. Once more factors (— l)^"~^f'-i trace Fermi or 
Bose character of Xaf). Using gjj(t) = ^'(-t)[.9fe (*)~5fe (01 
and utilizing p8)) and leads to final expression for the 
first correlation function in (|A1|) . Repeating considera- 
tion for the three other correlation functions in (jAip , and 
using the resulting expressions in (|XT|) leads to Eq. pT)) . 



APPENDIX B: GREEN FUNCTIONS IN THE 
LIOUVILLE SPACE 

Here we discuss properties of retarded and advanced 
Green functions in the Liouville space. We start from 
definitions ([50]) and ([5^ . Utilizing the property of the 
full unitary propagator 

<ci|e"*^*|^>=<cit|e-»£*|Bt -5>* (Bi) 
one can write 

< Xij\Ueff{t)\X,nn > =<C X ji\Uef f {t)\Xnm >* 

nrn 

\Ulff{t)\X,,-> (B2) 

where the second equality comes from definition of Her- 
mitian conjugate. Using (|B2|) in ([30|) one gets 

- ie{~t) « X,^^\Ueff{-t)\X,j » (B3) 

Note, that definitions ([50)1 and (|3ip lead to the usual 
Hcrmitian-type connection (jB3[) between retarded and 
advance Green functions = [Q'^]'^ . An alternative dcf- 
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inition 

5h,mnit) ^ -m « X,./K|e-*^*|l„„pf » (B4) 

mn 

g:^,^M = ^O^t) « XnrrJK\e'^'\X^Jp''^ » (B5) 
^ l6(-t) < X,,\Ulfj{-t)\Xran > 

would lead to Liouvillian conjugation^^ 5" — \Q^] ^ or 



APPENDIX C: EXPRESSION FOR C. 



Substituting (jCip and (|C2p into ([35)) and using standard 
Markov approximation 



Pab 



(ti) « ^ « a&|e*^»-^(*-*i)|cd » (C3) 



c,d 



one gets the (Markovian) Rcdficld quantum master equa- 
tion 



e/7 



We start from (|37p and p8|) and use free propagator 
in place of effective one. This leads to 



a'"'" (t) 



»j,m«\"/ -ie{t) <s: Xj^\e 

= -ie{t) < jle-'""' 



— iCAiti Y ^ 



-,(0) a 



{t) = i9{-t) <s: Xj,\e-'^^^''\X„ 



(CI) 

> 

(C2) 



id{-t) < j\e 



dt 



< ab\Cf,ff\cd > Pcd{t), 



where the generator for our model takes the form 



iC 



iC 



efff 



{a-b),{c-d) (b;a) ,(d;c) 



■EE 



_^jj(Wb + l)j^(Wb + l)^(Wi,)^(7Vi,)^> 



- ^JV„-l,JVe'5Ar,-l,Arrf(-l)^° X 



s^i ^ ri ^pj - ^(iV„-l,p;JV„,r),(7Vt-l,s<j;iVt,S6)V-^i -^j 



"•"f^ri ^SbJ ^Sdj ^ pj ^{N^-l,s,-Na,Sa),{Nb-l,p;Nb,r)\^i > 



2^\^ri ^si ^s^j ^pj ^{Ni,Sh;Nt + l,s),{Nh,p;Nh + l,r)\-^i > 

^{Nt-l,p;Nt,r),{Nt~l,s;Nt,St}\-^i ) 



P3 Saj ^{N^,p:Na + l,r),{Na.Sa;Na. + l,s)y^i I 



rt Set sj 



Pj '^{Na-l,s;Na,Sa),{Na-l,p;Na.r)\^-~ 



(iV„) _ ^.(A'a-l)- 



(C4) 
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where U^^^ are unitary transformations diagonalizing and are corresponding eigenvalues 

charge blocks H^^-* of the molecular Hamiltonian ([7]), 
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